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f^ ^ Abstract 

The present DAMA/LIBRA experiment and the former 

DAMA/Nal have cumulatively released so far the results ob- 

K> i tained with the data collected over 13 annual cycles (total exposure: 

^ • 1.17 ton X yr). They give a model independent evidence of the 

5t 1 presence of DM particles in the galactic halo on the basis of the 

DM annual modulation signature at 8.9 a C.L. for the cumulative 

exposure. 

1 Introduction 

The DAM A project is based on the development and use of low background 
scintillators, and several low background set-ups have been realized and used 
for various kinds of investigations [1]. In particular, the former DAMA/Nal 
and the present DAMA/LIBRA experiments at the Gran Sasso National 
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Laboratory of the INFN have the main aim to investigate the presence of 
Dark Matter particles in the galactic halo by exploiting the model indepen- 
dent DM annual modulation signature, originally suggested in the mid 80's 
in ref. [2]. In fact, as a consequence of its annual revolution around the 
Sun, which is moving in the Galaxy travelling with respect to the Local 
Standard of Rest towards the star Vega near the constellation of Hercules, 
the Earth should be crossed by a larger flux of Dark Matter particles around 
~ 2 June (when the Earth orbital velocity is summed to the one of the so- 
lar system with respect to the Galaxy) and by a smaller one around ~ 2 
December (when the two velocities are subtracted). It is worth noting that 
this signature has a different origin and peculiarities than the seasons on 
the Earth and than effects correlated with seasons (consider the expected 
value of the phase as well as the other requirements listed below). This 
annual modulation signature is very distinctive since the effect induced by 
DM particles must simultaneously satisfy all the following requirements: the 
rate must contain a component modulated according to a cosine function (1) 
with one year period (2) and a phase that peaks roughly around ~ 2 June 
(3); this modulation must only be found in a well-defined low energy range, 
where DM particle induced events can be present (4); it must apply only 
to those events in which just one detector of many actually "fires" {single- 
hit events), since the DM particle multi- interaction probability is negligible 
(5); the modulation amplitude in the region of maximal sensitivity must be 
<7% for usually adopted halo distributions (6), but it can be larger in case 
of some possible scenarios such as e.g. those in refs. [H Sj. This offers an 
efficient DM model independent signature, able to test a large interval of 
cross sections and of halo densities; moreover, the use of highly radiopure 
Nal(Tl) scintillators as target-detectors assures sensitivity to wide ranges of 
DM candidates, of interaction types and of astrophysical scenarios. 

It is worth noting that only systematic effects or side reactions able 
to simultaneously fulfil all the requirements given above (and no one has 
ever been suggested over more than a decade) and to account for the whole 
observed modulation amplitude might mimic this DM signature. 

The description, radiopurity and main features of the DAMA/LIBRA 
set-up are discussed in details in ref. [HIS]; moreover, this set-up has firstly 
been upgraded in September/October 2008 [7]. 

The cumulative DAMA/LIBRA exposure - after the new data release 
occurred at beginning of 2010 [7] - is 0.87 ton x yr (6 annual cycles), and 
cumulatively with DAMA/Nal the exposure is 1.17 ton x yr (13 annual 
cycles in total). 



2 The model independent result 

Several analyses on the model-independent investigation of the DM annual 
modulation signature have been performed in [7] as previously done in ref. 
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Figure 1: Experimental model- independent residual rate of the single-hit scintil- 
lation events, measured by DAMA/LIBRA-1,2,3,4,5,6 in the (2 - 6) keV energy 
interval as a function of the time ^6} 17j. The zero of the time scale is January 1** 
of the first year of data taking of the former DAMA/Nal experiment. The experi- 
mental points present the errors as vertical bars and the associated time bin width 
as horizontal bars. The superimposed curve is the cosinusoidal function behavior 
Acosuj{t — to) with a period T = — = 1 yr, with a phase to = f52.5 day (June 2"'') 
and with modulation amplitude, A, equal to the central value obtained by best fit 
over the whole data including also the exposure previously collected by the former 
DAMA/Nal experiment. The dashed vertical lines correspond to the maximum 
expected for the DM signal (June 2"''), while the dotted vertical lines correspond 
to the minimum. See refs. [6l [7] and refs. therein. 

[6] and refs. therein. In particular, Fig. [1] shows the time behaviour of 
the experimental residual rates for single-hit events in the (2-6) keV energy 
interval; as known, here and hereafter keV means keV electron equivalent. 
The hypothesis of absence of modulation in the data can be discarded [6l 
[7]. Moreover, when the period and the phase parameters as well as the 
modulation amplitude are kept free fitting the experimental residuals of 
Fig. [T] with the formula: Acoscj(t - to), values well compatible with the 
expectations for a signal in the DM annual modulation signature are found 
[6l[7]. In particular, the phase - whose better determination is obtained by 
using a maximum likelihood analysis [U [7] - is consistent with about June 
2"'^ within 2a. For completeness, we note that a slight energy dependence 
of the phase could be expected in case of possible contributions of non- 
thermalized DM components to the galactic halo, such as e.g. the SagDEG 
stream J8j and the caustics [9]. 

The data have also been investigated by a Fourier analysis, obtaining a 
clear peak corresponding to a period of 1 year; the same analysis in other 
energy region shows instead only aliasing peaks [6l[7]. 



The measured energy distribution has been investigated in other energy 
regions not of interest for Dark Matter, also verifying the absence of any 
significant background modulatioro- In particular, the measured rate in- 
tegrated above 90 keV, R90, as a function of the time has been analysed; 
fitting its time behaviour with phase and period as for DM particles, a mod- 
ulation amplitude compatible with zero is found excluding the presence of 
any background modulation in the whole energy spectrum at a level much 
lower than the effect found in the lowest energy region for the single-hit 
events [6l |7j . Similar result is obtained when comparing the single-hit resid- 
uals in the (2-6) keV with those in other energy intervals [SI [7] . It is worth 
noting that the obtained results already account for whatever kind of back- 
ground and, in addition, that no background process able to mimic the DM 
annual modulation signature (that is able to simultaneously satisfy all the 
peculiarities of the signature and to account for the measured modulation 
amplitude) is available (see also discussions e.g. in |6l[l0], refs. therein and 
later) . 

A further relevant investigation has been performed by applying the same 
hardware and software procedures, used to acquire and to analyse the single- 
hit residual rate, to the multiple-hit one. In fact, since the probability that 
a DM particle interacts in more than one detector is negligible, a DM signal 
can be present just in the single-hit residual rate. Thus, the comparison 
of the results of the single-hit events with those of the multiple-hit ones 
corresponds practically to compare between them the cases of DM particles 
beam-on and beam-off. This procedure also allows an additional test of 
the background behaviour in the same energy interval where the positive 
effect is observed. In particular, in Fig. [2] the residual rates of the single-hit 
events measured over the six DAMA/LIBRA annual cycles are reported, as 
collected in a single cycle, together with the residual rates of the multiple-hit 
events, in the considered energy interval. While, as already observed, a clear 
modulation, satisfying all the peculiarities of the DM annual modulation 
signature, is present in the single-hit events, the fitted modulation amplitude 
for the multiple-hit residual rate is well compatible with zero [7j. Thus, again 
evidence of annual modulation with proper features as required by the DM 
annual modulation signature is present in the single-hit residuals (events 



^In fact, the background in the lowest energy region is essentially due to "Compton" 
electrons, X-rays and/or Auger electrons, muon induced events, etc., which are strictly 
correlated with the events in the higher energy part of the spectrum. Thus, if a modulation 
detected in the lowest energy region would be due to a modulation of the background 
(rather than to a signal), an equal or larger modulation in the higher energy regions 
should be present. 



class to which the DM particle induced events belong), while it is absent in 
the multiple-hit residual rate (event class to which only background events 
belong). Since the same identical hardware and the same identical software 
procedures have been used to analyse the two classes of events, the obtained 
result offers an additional strong support for the presence of a DM particle 
component in the galactic halo. 
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Figure 2: Experimental residual rates over the six DAM A/LIBRA annual cycles 
for single-hit events (open circles) (class of events to which DM events belong) and 
for multiple-hit events (filled triangles) (class of events to which DM events do not 
belong). They have been obtained by considering for each class of events the data 
as collected in a single annual cycle and by using in both cases the same identical 
hardware and the same identical software procedures. The initial time of the figure 
is taken on August 7*''. The experimental points present the errors as vertical 
bars and the associated time bin width as horizontal bars. Analogous results were 
obtained for the DAMA/Nal data [TT]- See refs. [51[7]. 

The annual modulation present at low energy can also be shown by 
depicting - as a function of the energy - the modulation amplitude, Sm, 
obtained by maximum likelihood method over the data considering T =1 
yr and to = 152.5 day, as described in refs. [SlE]. In Fig. [3] the obtained 
Sm are shown in each considered energy bin (there AE = 0.5 keV). It can 
be inferred that positive signal is present in the (2-6) keV energy interval, 
while Sm values compatible with zero are present just above. In fact, the Sm 
values in the (6-20) keV energy interval have random fluctuations around 
zero with x^ equal to 27.5 for 28 degrees of freedom. All this confirms the 
previous analyses. 

The method also allows the extraction of the the Sm values for each 
detector, for each annual cycle and for each energy bin. Thus, following 
the procedure described in ref. p], we have also verified that the Sm are 
statistically well distributed in all the annual cycles and energy bins for 
each detector [6l [7j. Among further additional tests, the analysis of the 
modulation amplitudes as a function of the energy separately for the nine 
inner detectors and the remaining external ones has been carried out; the 
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Figure 3: Energy distribution of Sm for the total cumulative exposure 1.17 tonxyr. 
The energy bin is 0.5 keV. A clear modulation is present in the lowest energy region, 
while Sm values compatible with zero are present just above. In fact, the Sm values 
in the (6-20) keV energy interval have random fluctuations around zero with x^ 
equal to 27.5 for 28 degrees of freedom. See refs. [6l[7]. 



values are fully in agreement, showing that the effect is well shared between 
inner and external detectors [U [7] . 

Finally, releasing the assumption of a phase fo = 152.5 day in the max- 
imum likelihood procedure to evaluate the modulation amplitudes from 
the data of the 1.17 ton x yr exposure, one can alternatively write the 
signal as [6l \7\: Sik = S'o.fc + Sm,kCosu){ti - to) + Zm,ksinuj{ti - to) = 
So,k + ^m,fcCoscj(tj — t*). For signals induced by DM particles one would 
expect: i) Z^^k ~ (because of the orthogonality between the cosine and 
the sine functions); ii) Sm,k — Ym,k] iii) ** — ^o = 152.5 day. These condi- 
tions hold for most of the dark halo models; however, as mentioned above, 
slight differences can be expected in case of possible contributions from non- 
thermalized DM components, such as e.g. the SagDEG stream [8] and the 
caustics j;9j. 

Fig. m shows the 2a contours in the planes (5m, Zm) and (Ym,t*) for the 
(2-6) keV and (6-14) keV energy intervals. The best fit values are reported 
in [7]. Then, forcing to zero the contribution of the cosine function, the 
Zm values as function of the energy have also been determined by using the 
same procedure. The values of Zm are expected to be zero; by the fact, 
the x^ test supports such a hypothesis [7j. As in the previous analyses, 
an annual modulation effect is present in the single-hit events in the lower 
energy interval and the phase agrees with that expected for DM induced 
signals. These results confirm those achieved by other kinds of analyses. 

We further stress that sometimes naive statements were put forwards as 
the fact that in nature several phenomena may show some kind of period- 
icity. It is worth noting that the point is whether they might mimic the 
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Figure 4: 2cr contours in the plane {S„i, Zm) (left) and in the plane {Ym,t*) {right) 
for the (2-6) keV and (6-14) keV energy intervals. The contours have been obtained 
by the maximum likelihood method, considering the cumulative exposure of 1.17 
ton X yr. A modulation amplitude is present in the lower energy intervals and the 
phase agrees with that expected for DM induced signals. See refs. [SlIZ]- 



annual modulation signature in DAMA/LIBRA (and former DAMA/Nal), 
i.e. whether they might be not only quantitatively able to account for the 
observed modulation amplitude but also able to contemporaneously satisfy 
all the requirements of the DM annual modulation signature. The same is 
also for side reactions. This has already been deeply investigated e.g. in ref. 
[Bl El Ej and references therein. Some additional arguments have also been 
recently addressed in [inillS]. See also later. 

No modulation has been found in any possible source of systematics 
or side reactions for DAMA/LIBRA as well; moreover, no one is able to 
mimic the signature. Thus, cautious upper limits (90% C.L.) on the possible 
contributions to the DAMA/LIBRA measured modulation amplitude have 
been estimated and are summarized e.g. in Table of ref. [6]. 

Just as an example we recall here the case of muons, whose flux has 
been reported by the MACRO experiment to have an about 2% modulation 
with phase around mid- July [E]; recently, also LVD and Borexino results 
[T^l IT5] have been reported. We have already demonstrated that not only 
this effect would give rise in the DAMA set-ups to a quantitatively negligible 
contribution (see e.g. [H [71 [161 E] and refs. therein), but several of the six 
requirements necessary to mimic the DM annual modulation signature would 
also fail: e.g. muon would also induce modulation in the multiple-hit events 
and in the whole energy distribution, which is not observed. Moreover, even 
the pessimistic assumption of whatever hypothetical (even exotic) possible 



cosmogenic product - whose decay or de-excitation or whatever else might 
produce: i) only events at low energy; ii) only single-hit events; iii) no 
sizeable effect in the multiple-hits counting rate - cannot give rise to any 
side process able to mimic the investigated DM signature; in fact, not only 
this latter hypothetical process would be quantitatively negligible (see e.g. 
[HI E] and refs. therein), but in addition its phase would be (much) larger 
than July 15th, and therefore well different from the one measured by the 
DAMA experiments and expected by the DM annual modulation signature 
(~ June 2nd). In addition, the phase measured by DAMA experiments 
in each annual cycle is always around June 2nd, while the muon phase 
varies from year to year depending on the condition in the atmosphere. In 
particular, the value (146 it 7) days [7] measured by DAMA is 5.6 a distant 
from the LVD mean value (~ 185 days) and even more from the MACRO 
one; similar conclusions also hold for Borexino [15]. To be clear, let us 
note that even a phase value +3cr from the one measured for the DAMA 
observed effect - that is mid-June - cannot match the muon phase unless for 
one exceptional year; in fact, the maximum outer atmosphere temperature 
variation (and consequently the atmosphere density variation which causes 
the muon flux modulation) is typically not in that period at Gran Sasso 
location. In conclusion, any possible effect from muons is safely excluded on 
the basis of all the given quantitative facts (and just one of them is enough). 
There has been suggested in a recent paper |17) that environmental neu- 
trons (mainly thermal and/or epithermal) once captured by Iodine, might 
be responsible of the observed modulation through ^^^I decay, that pro- 
duces - among others - low energy X-rays/Auger electrons, provided any 
hypothetical modulation of the neutron impinging the sensitive part of the 
detectors inside the multi-component multi-ton shield. To be as clear as 
possible, we skip several comments about that paper and we just summa- 
rize few points. Firstly, we stress that - as already quantitatively discussed 
e.g. in |16t \TT\ El |7] - environmental neutrons cannot give any significant 
contribution to the annual modulation measured by the DAMA experiments 
lI and that ^^^I - if any - cannot mimic the DM annual modulation signature 
since some of its peculiarities would fail. Moreover, when the ^^^I decays in 
the EC channel, it produces low energy X-rays and Auger electrons, but - 



^For completeness, we also recall that the experimental set-ups, located deep under- 
ground, were equipped with a neutron shield made by Cd foils and polyethylene/paraffin 
moderator; moreover, a ~ 1 m concrete almost completely surrounds the installation 
acting as a further neutron moderator. The effectiveness of this shield has also been 
demonstrated e.g. in ref. [6] where a reduction larger than one order of magnitude has 
been measured for the thermal neutrons. 



since the ^^^I would be inside the Nal(Tl) detectors ~ the detectors would 
measure the total energy release of all the X-rays and Auger electrons emit- 
ted following the EC, that is the atomic binding energy either of the K shell 
(32 keV) or of the L shells (4.3 to 5 keV) of the ^^STe. The probability 
that so low-energy gamma's and electrons would escape a detector is very 
small; thus, we can conclude that: 1) the L-shell contribution would be a 
gaussian centered around 4.5 keV; but this is excluded by the DAMA data 
(see Fig. [3]). Moreover, the efficiencies to detect an event per one ^^^I decay 
are: 2 x lO^^, 6 x lO^^^ ^nd 2 x lO^^ in (2-4) keV, (4-6) keV and (6-8) 
keV respectively, as calculated by Montecarlo code. Thus, in addition, the 
contribution of ^^^I in the (2-4) keV - if any - would be similar to the one 
in the (6-8) keV, while the data exclude that; 2) the i^-shell contribution 
(around 30 keV) must be 8 times larger than that of L-shell, while no mod- 
ulation has been observed above 6 keV and, in particular, around 30 keV; 
3) the ^^^I also decays by /3~ with much larger branching ratio than EC 
and with (3~ end-point energy at 2 MeV. Again, no modulation has instead 
been observed in DAMA experiments at high energy [SI [7]. Moreover, the 
data collected by DAMA/LIBRA allow the determination of the possible 
presence of ^■^^I in the detectors. In fact, neutrons would generate ^^^I ho- 
mogeneously distributed in the Nal(Tl) detectors; therefore studying the 
characteristic radiation of ^^'^I decay and comparing it with the experimen- 
tal data, one can obtain the possible ^^^I concentration. The most sensitive 
way to perform such a measurement is to study the possible presence of 
the 32 keV peak (i^-shell contribution) in the region around 30 keV. This 
has already been done by DAMA - for other purposes - in ref. [18], where 
there is also no evidence of such a peak in the DAMA/LIBRA data; hence 
an upper limit on the area of a peak around 32 keV can be derived to be: 
0.074 cpd/kg (90% CL) p^. Considering the branching ratio of the process 
and the related efficiency, one can obtain a limit on possible activity of ^^^I: 
< 14.8;uBq/kg (90% CL). This upper limit allows to derive the maximum 
expected counting rate from ^^^I (see Fig. ^-Top)^ showing that in every 
case its contribution - if any - is negligible. The contribution is also neg- 
ligible even in the hypothetical case that the neutron flux had a 10% time 
modulation with the same phase, period as a DM signal; in fact, even in 
such a case the contribution to Sm is < 3 x 10"^ cpd/kg/keV at low en- 
ergy (see Fig. ^-Bottom) , that is < 2% of the DAMA observed modulation 
amplitudes. Therefore, for all the given arguments (and just one of them is 
enough), no role is played by ^^^I. 
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Figure 5: Top - Data points: cumulative low-energy distribution of the single-hit 
scintillation events measured by DAMA/LIBRA 5 above the 2 keV energy thresh- 
old of the experiment. Histogram (color online): maximum expected counting rate 
from ^^^I decays corresponding to the measured upper limit on ^^*I activity in the 
Nal(Tl) detectors: <14.8 ^Bq/kg (90% C.L.); see the data in ref. [18 and the text. 
Bottom - Data points: the DAMA measured modulation amplitude as a function 
of the energy. Histogram (color online) : maximum expected modulation amplitude 
multiplied by a factor 30 as a function of the energy from ^^*I decays correspond- 
ing to the measured upper limit on ^^*I activity given above and assuming the 
hypothetical case that the neutron flux had a 10% time modulation with the same 
phase, period as a DM signal. Therefore, the contribution from ^^^I is negligible 
and cannot mimic the Sm behaviour; in addition, it is worth noting that ^^^I never 
could mimic the DM annual modulation signature since some of its peculiarities 
would fail (e.g. ^^®I would induce modulation also in other energy regions, which 
is not observed). See text. 

3 Conclusions 



As regards the corollary investigation on the nature of the DM candidate 
particle(s) and related astrophysical, nuclear and particle physics scenar- 
ios, it has been shown that the obtained model independent evidence can 
be compatible with a wide set of possibilities as discussed e.g. in refs. 
[13 in El IS [201 [211 [221 [231 [M], in the Appendix of ref. [6] and in lit- 
erature (as e.g. [25], etc.); other possibilities are open. Moreover, as re- 
gards possible comparisons with direct, indirect and accelerators activities 



see e.g. discussions in [71 EJ [THl [26l [121 [27] , etc.. Here we just recall that no 
other experiment exists, whose result can be directly compared in a model- 
independent way with those by DAMA/Nal and DAMA/LIBRA. 

In conclusion, the six annual cycles of DAMA/LIBRA have further con- 
firmed the peculiar annual modulation of the single-hit events in the (2-6) 
keV energy region. The total exposure by the former DAMA/Nal and the 
present DAMA/LIBRA is 1.17 ton x yr and the confidence level for the 
observed effect is cumulatively about 9 a CL. The data satisfy all the many 
requirements of the signature and no systematics or side reactions able to 
mimic it (that is, able to account for the measured modulation amplitude 
and to simultaneously satisfy all the peculiarities of the signature) have been 
found or suggested by anyone over more than a decade. 

In near future new PMTs with higher quantum efficiency will be installed 
in order to lower the 2 keV energy threshold, increasing the experimental 
sensitivity and improving the corollary information on the nature of the 
DM candidate particle(s) and on the various related astrophysical, nuclear 
and particle physics scenarios. Moreover, it will also allow the investigation 
of other DM features, of second order effects and of several rare processes 
other than DM. A third generation R&D effort towards a possible Nal(Tl) 
ton set-up, DAMA proposed in 1996, has been funded by I. N.F.N, and is in 
progress. 
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